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A vital requirement for a quantum computer is the ability to locally address, with high fidelity, any of its 
qubits without affecting their neighbors. We propose an addressing method using composite sequences of 
laser pulses, which reduces dramatically the addressing error in a lattice of closely spaced atoms or ions, and 
at the same time significantly enhances the robustness of qubit manipulations. To this end, we design novel 
high-fidelity composite pulses for the most important singlc-qubit operations. In principle, this method allows 
one to beat the diffraction limit, for only atoms situated in a small spatial region around the center of the laser 
beam are excited, well within the laser beam waist. © 2013 Optical Society of America 



Scalable quantum computers depend critically on the 
ability to perform local addressing of their individual 
qubits [1]. In a Paul ion trap, which is one of the most 
promising scalable platforms for the future quantum 
computer [2], local addressing is the ability to operate 
on a single ion by focused laser light while keeping the 
neighboring ions unaffected. When the number of ions 
increases, the distance between them diminishes and lo- 
cal addressing becomes one of the principal experimental 
challenges. For example, in a recent experimental demon- 
stration of the Toffoli gate [3] most of the error was at- 
tributed to addressing error, as the neighboring ions were 
seeing 7% of the central Rabi frequency. 

In this Letter, we propose a method for high-fidelity 
local addressing applicable to various types of atomic 
qubits: trapped ions, atoms in optical lattices, quantum 
dots, etc. To this end, we present new narrowband (NB) 
and passband (PB) composite pulses 1 , which are spe- 
cially designed for local addressing. The excitation pro- 
files of such pulses allow one to manipulate only a sin- 
gle qubit, as the outer parts of the spatial laser beam 
profile practically do not excite its neighbors, although 
the latter may be subjected to significant laser intensity. 
Moreover, with a PB pulse one enhances the robustness 
of qubit manipulations, thereby eliminating errors due 
to imperfectly calibrated and fluctuating laser intensity 
and laser beam pointing instability. 

The technique of composite pulses was introduced in 
nuclear magnetic resonance (NMR) [4-6] as a powerful 
tool for manipulation of spins by magnetic fields. A com- 
posite pulse compensates the imperfections of a single 
pulse, which is the traditional tool used to drive a quan- 
tum transition, and it consists of a sequence of pulses, 
each with a well-defined phase. The composite phases are 
determined from the conditions imposed on the desired 
overall excitation profile. In particular, in a NB compos- 
ite pulse only the qubits seeing pulse areas within a nar- 
row range around some value A are subjected to trans- 



formation, while qubits seeing areas outside this range 
remain unaffected in the end of the composite sequence. 

Most known composite pulses, however, are inappro- 
priate for local addressing due to the sidebands and the 
slowly vanishing tails in their excitation profiles. The few 
composite pulses that allow local addressing [7] are very 
sensitive to variations in the pulse area. We present here 
a systematic method, which allows us to construct high- 
fidelity composite pulses, which eliminate these draw- 
backs and which can be made robust and accurate to 
any desired order. 

A two-state quantum system ip 1 -h> -02, subjected to 
an external coherent electromagnetic field, obeys the 
Schrodinger equation, \hdtc(t) = H(t)c(t). Here the 
vector c(t) = [ci(£), C2(t)] T contains the two prob- 
ability amplitudes and the Hamiltonian is H(i) = 
(h/2)n(t)c~ iD ^\4> 1 )(i J2 \+h.c., with D(t) = £ A(t')di', 
where A = luq — uj is the detuning between the laser car- 
rier frequency u and the Bohr transition frequency ojq, 
and Q(t) is the Rabi frequency. The amplitudes at the 
end of the interaction c(t{) are obtained from the initial 
ones with the propagator U: c(tf) = U(if,ti)c(ii), which 
can be parameterized with the complex Cayley-Klein pa- 
rameters a and b (obeying |a| 2 + \b\ 2 = 1), 



U = 



a 

-b* 



(1) 



For resonance (A = 0), the Schrodinger equation has an 
exact solution regardless of the shape of O(i) and the 
parameters a and b are determined only by the pulse 

= -isin(A/2). The 
sin 2 (A/2). A con- 
stant phase shift </> in the driving field, fi(t) — > f2(t) e 1 ^, 
is mapped onto the propagator as 



area A = f* f fl(t)dt: a = cos(A/2), b 
transition probability is p = \b\ 2 = 
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1 We follow the usual NMR terminology, which is related to the 
features of the excitation profile, rather than the radiation. 



We assume for simplicity, and possible experimental 
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Narrowband sequences 


TV J ( A \ 

N 2n +i{A) 


A 


1 ( ± ± ± ± \ 

phases [fa; <p 3 ; </> 4 ; . . . ; <p n +i) 


N 5 (tt) 


TT 


(0.839; 1.420) 


N 9 (tt) 


TT 


(0.426; 1.490; 0.858; 1.300) 


Ni 3 (tt) 


TT 


(1.103; 0.876; 0.154; 1.708; 






1.020; 0.229) 


lN 2 i(,7rJ 


7T 


(i.U(o; u.yiy; u.ioi, i.ooi, 






1 156- 721- 096- 1 521- 






0.812; 1.954) 


N 7 (f) 


3tt/7 


(0.471; 1.196; 1.315) 




tt/2 


(0.577; 1.161; 1.573) 




3tt/8 


(1.532; 0.800; 0.698) 


N 7 (\/27t) 


3tt/4 


(1.505; 0.823; 0.609) 


Passband sequences 


P2n+l{A) 


A 


phases (<p 2 ; 03; 04 5 • ■ • ; <t>n+i) 


P 7 (7T) 


7T 


(0.508; 1.337; 1.083) 


Pl7(7T) 


TT 


(1.235; 0.721; 0.934; 0.126; 






1.872; 1.515; 0.873; 0.217) 


P 9 (f) 


3tt/5 


(1.270; 1.106; 0.464; 0.053) 


Pir(f) 


2tt/3 


(0.459; 0.097; 0.302; 1.445; 






0.829; 1.324; 1.290; 0.995) 


P »W 


3tt/5 


(0.676; 0.87; 1.503; 1.836) 




3tt/5 


(0.356; 1.517; 0.957; 1.023) 


P 9 (\/2tt) 


3tt/5 


(1.909; 1.197; 0.861; 0.660) 



Table 1. Phases <fik (in units tt) for some NB (Njv) and 
PB (Pjv) sequences of N = 2n + 1 phased resonant 
pulses of area A: A A < / >2 A < f >3 ■ ■ ■ A<f )ri+1 ■ ■ ■ A^ 3 A^, 2 A . We 
set rii = n , n 2 = for all NB pulses iV 2n+ i(„4); ri\ = 2, 
n 2 = 1 for P 7 (-4) and P 9 (-4); n\ = 6, n 2 = 3 for P^tt); 
n\ = 4, 7i 2 = 3 for Pn(n/2). For all sequences we 
set ?i3 = 0. For A ^ tt we also impose Eq. (4a). The 
pulse area is A = tt for N 5 (ir), N 9 (ir), Ni 3 (it), N 2 i(tt), 
P 7 {n) and Pi 7 (tt); A = 3tt/7, tt/2, 3tt/8 and 3tt/4, re- 
spectively, for N 7 (n/2), N 7 (n/V2), N 7 {n/{2-f2)) and 
N 7 (V2n); A = 2tt/3 for Pi 7 (tt/2); A = 3tt/5 for P 9 (7r/2), 
P 9 (7r/V2), Pg(7r/(2V2)) and P 9 (V2tt). 

convenience, that all pulse areas are equal 2 . A sequence 
of N pulses A r f >k , each with area A and phase 4>k , pro- 
duces the propagator 



uW=U(4 fe )U(4 fa _,)-U(^). 



(3) 



We consider here an odd number of pulses, N = 2n + 1, 
although this restriction is not essential. We consider 
composite sequences, which are symmetric with respect 
to reversal of pulses, i.e. the phases should obey 4>k = 
4>N+\-k (k = 1,2, ...,n); this "anagram" condition an- 
nuls the imaginary part of the propagator element . 
Since the overall phase of the sequence is irrelevant, but 
only the relative phases of the pulses matter for the dy- 
namics, we set 6\ = 4>n = 0; hence we are left with n in- 



2 Allowing for different pulse areas, and thereby seemingly for 
more free parameters, does not seem to be advantageous. We have 
found through extensive examination that the total pulse area (and 
the pulse duration) of the composite sequence is not reduced. 
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Fig. 1. Top frames: Excitation probability p in a laser 
field with a Gaussian spatial profile (grey-shaded) with 
FWHM of Rabi frequency £, vs the distance from its 
center for various composite pulses from Table 1: N5(7r), 
P 7 (tt), Nai(Tr) and P 17 (tt) (left) and N 7 (tt/2), P 9 (tt/2) 
and Pi 7 (7t/2) (right). The excitation profiles of single tt 
and 7t/2 pulses are shown too. Bottom frames: Deviation 
Po — p, with po = 1 (left) and po = 0.5 (right) being the 
desired excitation probability. 

dependent phases, which are treated as free parameters. 
For NB pulses we require transition probability p m for 
pulse areas in the vicinity of A = (a flat bottom); for 
PB pulses we also require p w const around the desired 
overall area A (a flat plateau). Thus, the phases for a 
NB sequence are derived from the conditions 
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A=A = COs(A/2), 
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,2m), 



(4a) 
(4b) 



with m = n—1 and d\ = d k /dA k . For a n 2 -fold flat-top 
PB sequence we add the conditions 



[d k A Uip] A=A = (k = 1,2,.. 



,n 2 ) 



(5) 



with m +n 2 + 1 = n. For a sequence with a target phase 
tp we must also have 

[argC/ 2 i U=A = V- (6) 
This phase can be stabilized with PB pulses, for which 

[axgd A u£ ) ] A = A = (fc = l,2,...,n 3 ). (7) 

Thus a phase-stable PB sequence consists of n — ni + 
n 2 + ri3 + 2 pulses. For target area A = tt, Eqs. (4a) and 
(5) for even k and (7) for odd k are fulfilled identically. 

Using the method described above, we have con- 
structed new NB and PB composite pulses; examples 
are presented in Table 1 for the most often used areas. 
Figure 1 illustrates the excitation profiles for different 
composite pulses of target area tt (left frames) and 7r/2 
(right frames). At the center of the laser spot (origin) the 
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Sequences 


A 


phases (&>>; <f>3\ 4>i\ ■ ■ ■ 5 4>n+i) 


N 7 (tt, f ) 
Pii(T,f ) 

N 9 (f,i) 

Pl3(f , |) 

P / 7T 37T \ 

F 13(^5, -J" ) 


7T 
7T 

3tt/4 
3tt/4 

3tt/5 
3tt/5 


(1.256; 0.792; 0.072) 
(0.221; 1.109; 0.753; 1.304; 
1.878) 

(1.074; 0.935; 0.173; 1.562) 
(0.959; 1.048; 0.367; 1.967; 

1.511; 0.860) 
(1.326; 0.958; 0.137; 0.791) 
(0.183; 0.978; 1.421; 0.769; 

1.924; 1.916) 



Table 2. Phases <pk (hi units ir) for some NB (Njv(-A, ip)) 
and PB (Pn(A 7 p)) sequences of N = 2n+l phased reso- 
nant pulses of area A: AgA^A^ ■ ■ ■ A<j, n+1 ■ ■ ■ A^A^Aq, 
which produce phased rotations of angle p at area A. 
We set n 2 = n 3 = for all NB pulses N 2n +i{A, ip). We 
have rii = 1, n 3 = 4 for Pn(A,ip); n 2 = n 3 = 1 for 
P13 (A, (p) . For all sequences we set ri\ = 2 and we im- 
pose Eq. (6). For A 7^ ir we also impose Eq. (4a). The 
pulse area is A ~ ir for Nt(tt,3tt/2) and Pu(tt, 3ir/2); 
A = 3tt/4 for N g (n/2, tt/2) and Pi 3 (tt/2, tt/2); A = 3tt/5 
for N 9 (w/V2,3w/2) and P 13 (ir/V2, 3tt/2). 



probability of excitation is p = sin 2 (.4/2), whereas at 
the wings it naturally decreases. The logarithmic scale 
allows us to examine the fidelity of the profile against 
the 10 -4 quantum information benchmark [1]. We as- 
sume that the laser beam is a spot with a full-width- 
at-half-maximum (FWHM) of intensity £/V2; this im- 
plies FWHM of Rabi frequency £ for single-photon tran- 
sitions. Remarkably, suppression of unwanted neighbor 
excitation below the 10" 4 benchmark is achieved at up 
to 15% of the peak Rabi frequency (corresponding to 
distance 0.83£ from the center of the spot) with the N5 
pulse, 27% (distance 0.70£) with P 17 , and 48% (distance 
0.51£) with N21. In this manner we can beat the diffrac- 
tion limit as the excitation is localized in a spatial range, 
which, with sufficiently many ingredient pulses, can be 
made smaller than the beam waist. The physical reason 
for this suppression is the destructive interference of the 
ingredient pulses in the composite sequence. 

The bottom frames of Fig. 1 show the infidelity of the 
target qubit itself and reveal that PB composite pulses 
can greatly enhance the robustness of manipulation of 
the target qubit without losing selectivity. For target 
area A = it (left frame) an infidelity of 10~ 4 is encoun- 
tered at offset 0.05£ for a single pulse, while the admis- 
sible offset reaches 0.18£ for P 7 (tt) and 0.21£ for Pi 7 (tt). 

Our scheme can be further exploited, in conjunction 
with Eqs. (6) and (7), to generate various NB and PB 
sequences of specified target phases ip at the expense 
of additional pulses. Table 2 lists a set of composite se- 
quences, which produce experimentally relevant phased 
rotations. Larger sequences achieve stabilization of the 
phase p. Figure 2 illustrates the robustness of the com- 
posite phase ip = Sir/ 2 [3] produced by two of our phase- 
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Fig. 2. Phase deviation \<f> — 3tt/2| vs pulse area devi- 
ation of the composite pulses Nj(tt, ^-) and Pu(ir, 4^), 
shown in Table 2. The dashed curves are for frequency 
mismatch A = 0.001/T. The horizontal line is the 
threshold above which the infidelity exceeds 10 -4 . 

stabilized sequences from Table 2. Remarkably, one can 
perform high-fidelity local addressing even with pulse 
area deviation of 20%. In addition, we have found that 
pulse area noise with relative amplitude 5% introduces 
an absolute error in the target phase of 2.5 x 10 _3 tt. 

To conclude, the new composite sequences designed for 
high-fidelity local addressing in a lattice of closely spaced 
qubits, are of potential application to Paul ion traps [7,8] 
and ultracold atoms in optical lattices [9] . They allow one 
to stabilize both the rotation angle and the phase of the 
desired qubit rotation. This technique can be adapted 
to addressing on the vibrational sidebands, which should 
allow to construct high-fidelity two-qubit operations. 

This work has been supported by the European Com- 
mission project FASTQUAST and the Bulgarian NSF 
grants VU-I-301/07, D002-90/08 and DMU02-19/09. 
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